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Abstract 

This paper presents a variable neighborhood search (VNS) algorithm for solv¬ 
ing bandwidth coloring problem (BCP) and bandwidth multicoloring problem 
(BMCP). BCP and BMCP are generalizations of the well known vertex color¬ 
ing problem and they are of a great interest from both theoretical and practical 
points of view. Presented VNS combines a shaking procedure which perturbs 
the colors for an increasing number of vertices and a specific variable neighbor¬ 
hood descent (VND) procedure, based on the specially designed arrangement of 
the vertices which are the subject of re-coloring. By this approach, local search 
is split in a series of disjoint procedures, enabling better choice of the vertices 
which are addressed to re-color. The experiments show that proposed method 
is highly competitive with the state-of-the-art algorithms and improves 2 out of 
33 previous best known solutions for BMCP. 

Keywords: Bandwidth Coloring, Bandwidth MultiColoring, Frequency 
Assignment, Variable Neighborhood Search, Variable Neighborhood Descent 


1. Introduction 

Vertex coloring problem (VCP) and its generalizations belong to a well 
known and widely researched class of graph coloring problems and therefore 
represents a major challenge for most researchers in the field of combinatorial 
optimization. Various generalizations and variants of the VCP have been re¬ 
searched over the years and there are thousands of scientific papers proposing 
various methods for solving VCP and its generalizations. 

In VCP, one needs to color the vertices of the graph in such a way that 
adjacent vertices must be colored with different colors and the aim is to mini¬ 
mize the number of used colors. During the years, as one of the most studied 


‘Corresponding author. phone/fax:+38751319142 
Email addresses: matic.draganSginail.com (Dragan Matic), jkratica8gmail.com (Jozef 
Kratica) 


Preprint submitted to Els 


May 25, 2015 





NP-hard combinatorial optimization problems, VCP has undergone many gen¬ 
eralizations. This paper deals with two generalizations of the VCP: bandwidth 
coloring problem (BCP) and bandwidth multicoloring problem (BMCP). 

BCP is a straightforward generalization of the VCP, where for each two 
adjacent vertices u and v the distance, d{u,v) is imposed and the difference 
between two colors assigned to u and v must be larger than or equal to the 
defined distance. The task is to color the vertices with the smallest number of 
colors. Formally, for a graph G = {V, E) with positive integer distance function 
d{u,v), {u,v) e E, the objective of the BCP is to find the coloring c such that 
for each pair of adjacent vertices u and v, \c{u) — c(v) | > d(u, v) and the number 
of colors is minimized. Obviously, if the distance between any pair of adjacent 
vertices is equal to 1, the BCP is brought down to the VCP. 

BMCP generalizes BCP by including the multicoloring of the vertices. For 
each vertex v in the input graph a positive integer weight w{v) is introduced, 
giving the requirement how many colors must be assigned to that vertex. Addi¬ 
tionally, distance between the vertex to itself is also given, holding the condition 
inherited from BCP. Formally, for a given graph G = (V, E) with positive inte¬ 
ger distance function d{u, v), (u, v) G E, with the condition that u can be equal 
to V (loops are allowed) and for the given positive integer value w{u), for each 
vertex u G V, the aim of the BMCP is to find the coloring c such that each 
vertex u is colored with w{u) colors, for each pair of adjacent vertices u and w, 
|c(m) — c(?;)| > d{u,v) (including loops) and the number of colors is minimized. 
In a case when w{u) = 1 for each vertex u, BMCP is reduced to the BCP. 

BMCP can be converted into the BCP, by splitting each vertex v into a 
clique of cardinality w{v). Each edge in the clique is assigned the distance 
d{u,u), corresponding to the distance of the loop edge of the vertex u in the 
original graph. By this approach, each instance of BMCP, with n vertices, is 
transformed to the instance of the BCP, having w{i) vertices. This fact 
leads to the approach of constructing the algorithm for solving only the BCP 
which, after the explicit or implicit construction of the appropriate graph, can 
also be applied to solve BMCP. Both BCP and BMCP are NP-hard problems, 
since they generalize VCP, which is NP-hard in a general case. 

Like many other graph - based problems, VCP and its generalizations enjoys 
many practical applications. 

For example, it is well known that timetabling problems can be interpreted 
as graph coloring problems. Timetabling problems typically includes the task 
of assigning timeslots to the events. In university timetabling problems, events 
(lectures or exams) are interpreted as vertices, constraints by edges and timeslots 
by colors. According to [l|, whole class of university timetabling problems can 
be roughly divided into the subclasses: course timetabling problems (CTP) and 
examination timetabling problems (ETP). For example, in CTP, lessons are 
vertices and two vertices are adjacent if corresponding lessons are taught by the 
same professor. The task is to assign the timeslots to the lessons (i.e. assign 
the colors to the vertices), in such a way that two lessons given by the same 
professor must not be taught at the same time. In ETP, exams are interpreted 
by the vertices and two vertices are adjacent if there is at least one student 
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passing the both exams. If a new constraint is introduced by including the 
required time distance between two exams, the coloring problem which arises 
from this approach becomes BCP. Additional constraints can be placed on the 
each vertex, for example, if one exam has to be organized for many groups of 
students. In that case, the corresponding coloring problem is the BMCP. 

Recent approaches for solving timetabling problems, based on using graph 
coloring techniques, involves the hybridization of the methods automated by 
the hyper-heuristic methods: a generic hyper-heuristic approach of various con¬ 
structive heuristics for solving the VCP is presented in Q. In a consecutive 
work 0, the authors presented a random iterative graph based hyper-heuristic, 
which adaptively hybridise two graph colouring heuristics at different stages of 
solution. Another recent hyper-heuristic approach utilizes the hierarchical 
hybridizations of four low level graph coloring heuristics: largest degree, sat¬ 
uration degree, largest colored degree and largest enrollment. A constructive 
heuristic for finding a feasible timetable is recently presented in Q . Much more 
information and a recent general survey of graph coloring algorithms can be 
found in Q. 

An intensive grow of the usage of the wireless communication protocols, like 
in mobile telephony, satellite communication, wireless sensor networks, radio 
and TV broadcasting etc. leads to the development of various strategies ap¬ 
plied to organize and manage the frequencies of the signals: the handling of 
interference among radio signals, determining the availability of frequencies, as 
well as including various optimization criteria for achieving better performances 
of the overall system. The optimization problems which arises from the task of 
managing frequencies fall in the class of so called frequency assignment problems 
(FAP) and according to the special needs of a concrete problem, many variants 
of the FAP appears. A common feature of most of the variants includes the dis¬ 
tance constraint imposed on pairs of frequencies, in order to avoid or reduce the 
interference between close communication devices. In other words, two commu¬ 
nication points which are close enough (as so can interference each other) must 
be assigned with the enough different frequencies. Without any additional con¬ 
straint, this frequency assignment problem is simply called feasible frequency 
assignment problem (F-FAP) and it is obvious that this F-FAP corresponds to 
the finding of feasible coloring of the corresponding graph. 

In connection with FAP problems. Hale Q introduced so called T-coloring 
of the graph. For a given graph G = (V, E) and for each edge {u, w} £ E, 
we associate a set Tu,v of nonnegative integers, containing 0. A T-coloring of 
G is a function (an assignment) c : V{G) —>• N of colors to each vertex of 
G, so that if u and v are adjacent in G, then |c(u) — c(v)| is not in Tij. In 
simple words, the distance between two colors of adjacent vertices must not 
belong to the associated set T. The span of a T-coloring c is the difference 
between the smallest and the highest color in c. The optimization version of T- 
coloring is finding the minimum span of all possible T - colorings. It is obvious if 
Tu,v = {0} for each edge {u, u}, then problem becomes VCP. A generalization of 
T-coloring is a set T-coloring problem Q, which introduces the multiple coloring 
of vertices in the graph. For each vertex v, the set T-coloring problem includes 
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the assignment of a nonnegative integer 5^, providing the information how many 
colors need to be assigned to v. Also, for each vertex v a set is introduced. 
The constraints for each pair of adjacent vertices, previously introduced in T- 
coloring problem are extended with the constraints related to each vertex v: 
the distance of any two numbers (colors) assigned to v must not be in Ty^y. 
BCP is a restriction of T-coloring, where the constraint on adjacent vertices is 
replaced by the proper condition of BCP: |c(m) — c(u)| > t{u,v), where t(u,v) 
is a numerical value. Similarly, BMCP can be considered as an instance of the 
set T-coloring problem, where the set Ty^y is replaced by the numerical value, 
corresponding to d(y,v). A survey of the results and problems concerning with 
T-colorings can be found in 

According to the other special needs of a concrete problem, various subprob¬ 
lems of FAP can arise [l^ : 

• The Maximum Service (Max-FAP) and Minimum Blocking Frequency As¬ 
signment Problems (MB-FAP) : if feasible solutions to the F-FAP are not 
available or are difficult to find, a partial solution that assigns as many 
frequencies as possible to the vertices can be considered. This problem is 
similar with registry allocation problem. 

• The Minimum Order FAP (MO-FAP): Frequencies are assigned in such a 
way that no unacceptable interference occurs, and the number of different 
used frequencies is minimized. This problem is a direct generalization of 
the graph coloring problem. 

• The Minimum Interference Frequency Assignment Problem (MI-FAP): fre¬ 
quencies from a limited number of available frequencies are assigned in 
such a way that the total sum of weighted interference is minimized; in 
fact in this problem, we minimize the sum of the penalties incurred by the 
frequency choices 

• The Minimum Span Frequency Assignment Problem (MS-FAP): the prob¬ 
lem is to assign frequencies in such a way that the interference between 
the points is avoided, and the difference between the maximum and min¬ 
imum used frequency, the span, is minimized. This problem is equivalent 
to the BCP. In a case when the feature that one point is to be assigned 
with more than one frequency, the problem becomes the BMCP. 

• Other variants with some additional constraints 


More applications, as well as other discussions considering graph coloring 
and its generalizations is out of the paper’s scope and can be found for example 

ilQ- 


The rest of this paper is organized as follows. The next section recapitu¬ 
lates previous work regarding BCP and BMCP. VNS approach is described in 
details in the section [3l Section |4] contains experimental results obtained on the 
instances from the literature, while the last section contains conclusions and 
future work. 
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2. Previous work on solving BCP and BMCP 


As stated above, BCP and BMCP, as well as MS-FAP has been intensively 
solved by a large number of successful methods. In 2002, in order to encour¬ 
age research on computational methods for solving graph coloring problems, a 
’’Computational Symposium on graph coloring and its generalization” was or¬ 
ganized. The Symposium included the following topics: exact algorithms and 
heuristic approaches for solving the graph coloring problems, applications and 
instance generation, as well as methods for algorithm comparison [l5 |. In or¬ 
der to make a fair comparison of the proposed methods, for testing the solving 
the bandwidth (edge weights) multicoloring (node weights) or both graph col¬ 
oring problems, a standard set of test instances was suggested. This set is also 
used in this paper. At the Symposium, several successful heuristic methods 
were proposed. A problem-independent heuristic implementation called Dis- 
cropt, designed for ’’black box optimization”, was adapted to graph coloring 
problems and provided a good test of the flexibility of the system 1^. At 


the same symposium, Prestwitch ( 0 ) proposed a heuristic based on the hy¬ 
bridization of local search and constraint propagation. In the consecutive con¬ 
tribution [16|, the author extended his previous work by adding the constraint 
programming technique of forward checking in order to prune the colouration 
neighbourhoods. Hybrid methods using a squeaky-wheel optimization, com- 
bined with hill-climbing and with tabu heuristic, were described in 17| and 
[isj l . Chiarandini et al. in [l^ presented an experimental study of local search 
algorithms for solving general and large size instances of the set T -colouring 
problem. 

Malaguti and Toth (0) proposed a successful method for solving BCP and 
BMCP, which combines an evolutionary algorithm with a tabu search procedure. 
Like the method used in [0, Malaguti and Toth’s algorithm starts with the 
construction of an initial solution with a greedy approach. After that, it tries 
to improve the starting solution by reducing by one unit the number of colors 
used. Marti et al. ^ proposed the memory-based and memory-less methods 

oroblem, based on tabu search and GRASP, 
uses Multistart iterated tabu search (MITS) 


to solve the bandwidth colouriM 
Paper presented by Lai and Lii [22[ 
algorithm, which integrates an iterated tabu search with multistart method 
and a problem specific operator designed for the perturbation. This method 
uses previous known best results as a starting result and tries to improve it by 
decreasing it by one, repeating this process until no legal coloring can be found. 
The tabu search proposed in 22| is successfully combined with a path relinking 
algorithm (paper available at ArXiv 0). Computational results demonstrate 
that the proposed algorithm outperforms previous methods, improving 15 out of 
66 instances and matches 47 previous best known results. Lastly, a very recent 
work presented by Jin and Hao 2^ uses a learning-based hybrid search for 
solving BCP and BMCP. The algorithm firstly builds feasible partial colorings, 
which are further improved in a local search phase. In the construction phase, a 
learning-based guiding is used to determine the next vertex for color assignment, 
while in the local search phase tabu search technique is used to repair the 


5 








solution. This approach further improves the best solutions for 14 common 
used instances. 

From the last paragraph it is evident that in recent years the BCP and 
BMCP have been intensively solved by many highly efficient heuristics. As a 
consequence, the solutions for the benchmark data got by these methods are of 
a very high quality. In the next section, we describe the variable neighborhood 
search, which can further improve two of these intensively solved benchmark 
instances. 


3. VNS for solving BCP and BMCP 

This section presents the VNS for solving BCP and BMCP. Recall that each 
BMCP instance is implicitly transformed to an instance of BCP, by replacing 
each vertex v of the weight w(v), by the clique of the size w(v). Therefore, only 
the algorithm for solving BCP is described, but it is also applied on the BMCP, 
after the mentioned transformation of the instances. 

Variable Neighborhood Search (VNS) algorithm was originally described by 
Mladenovic and Hansen ([1^, [l^). In recent years, VNS has been proven as 
a very effective and adoptable metaheuristic, used for solving a wide range of 
complex optimization problems. The basic strategy of the VNS is to focus the 
investigation of the solutions which belong to some neighborhood of the current 
best one. In order to avoid being trapped in local suboptimal solutions, VNS 
changes the neighborhoods, directing the search in the promising and unexplored 
areas. By this systematic change of neighborhoods, VNS iteratively examines a 
sequence of neighbors of the current best solution, following the approach that 
multiple local optima are often in a kind of correlation, holding the ’good parts’ 
of the current best solution and trying to improve the rest of it. 

Many successful implementations of the standard VNS, as well as its many 
variants, prove that this successive investigation of the quality of the current 
solution’s neighbors can lead to better overall solutions. Standard VNS usually 
imposes two main procedures: shaking and local search (LS). Shaking procedure 
manages the overall system of the neighborhoods and in each iteration suggests 
a new point (potential solution) from the current neighborhood. In order to 
better widen the search, shaking procedure often uses the neighborhoods of the 
different cardinality. More precisely, for the given numbers kmin and kmax, 
a system of neighborhoods ..., is constructed. For each 

value k G [kmin, kmax], a (usually random) solution from the neighborhood Nk 
is chosen, which is the subject of further possible improvement inside the LS 
procedure. 

LS is trying to improve the suggested solution, by investigating the other 
solutions in its neighborhood, usually formed by some minor changes of it. 
Local search is usually implemented by using either best improving strategy 
or first improving strategy. While in the best improving strategy local search 
investigates all neighbors, keeping the current best one as a new solution, the 
first improving strategy stops when the first improving neighbor is found. 
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In the proposed VNS, we enhance the basic LS approach by ’’splitting” 
the local search in a disjoint union of smaller ones, enabling the algorithm to 
’’switch” between different neighborhoods inside the LS procedure. The current 
neighborhood is analysed until the solution is improved. If so, it is restarting 
from the first neighborhood, otherwise, goes into the next one. This approach 
requires the introduction of a specific variable neighborhood descent (VND) 
procedure which manages the examination of the neighborhoods, allowing the 
algorithm to firstly analyse the best suitable one. The description of the VND 
is in details explained in the section [3.41 

The optimization process of the VNS algorithm finishes when the stopping 
criterion is achieved, usually given by the maximum number of iterations or 
maximum time allowed for the execution and the latter is the case in this paper. 

Some recent successful implementation of the VNS using the VND approach 
can be found for example in . 

The overall VNS algorithm is shown on the figure [TJ The algorithm inputs 
the following data: distance constraint matrix, values kmin and k^ax, denoting 
the minimal and maximal neighborhood structures, timcmax - maximal allowed 
execution time and the value Pmove, representing the probability of shifting 
from one solution to another, in a case of equal objective functions. After the 
data input, VNS starts with a greedy heuristic (described in the subsection l3.1l) . 
which gives the upper bound (UB) for the total number of colors. After the UB 
is determined, the initial value of k* is set to UB (legal coloring with UB colors), 
and the starting solution of the VNS is constructed by the procedure Init(). 
Initial solution uses one color less than the value UB. This procedure, together 
with the objective function, is described in the subsection l3.2l The minimization 
process is performed in the shaking procedure (subsection 13.31) and the VND 
procedure, which is, together with the Compare procedure described in the 
subsection l3.4l During the minimization process inside the VND procedure the 
algorithm is trying to improve the solution to the feasible one. If this situation 
happens, the value fc* is decreased by one, that legal coloring is remembered 
and the algorithm repeats the search process by decreasing the total number of 
colors. The algorithm stops when the maximum execution time is reached. The 
result of the algorithm is the value k*, i.e. the fc*-coloring of the given graph. 


3.1. Constructive heuristics 

Before the minimization process is started, an initial solution (legal fc-coloring, 
for some k) needs to be constructed. Our algorithm uses the gr eedy approach 
similar to the greedy algorithm used and minutely described in [l8l [SOj . This 
greedy algorithm takes a sequence of ’split nodes’, greedily assigning colors to 
them. For each node a set of ’forbidden’ colors is firstly formed and after that 
the algorithm chooses the smallest color not belonging to that set. 

In the literature, other approaches for getting starting solutions are also 
used. For example, Malaguti and Toth in [20|, considered several greedy algo¬ 
rithms proposed for solving VCP: sequential greedy algorithm (SEQ), as well as 
another greedy approach - DSATUR from 3l|, similar to SEQ, but one which 
dynamically chooses the vertex to color next, i.e. the one which minimizes a 
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Data: Graph, Distance constraint matrix, kmin, k^ax, timemax^ Pmove 
Result: number k* 
xs ^ GreedyHeuristic(); 

X ^RemoveLastGolor(a;s); 

k kmin , 

while time < timemax do 
x' = Shaking (a;, fc); 

a;"=VND(a;,a;'); 

if Comparefx”,x,Pmove) then 

I X ■«— x"; 

else 

if fc < kmax then 
I k i — A: -f 1 5 
else 

I k i kmin j 

end 

end 

end 


Figure 1: VNS pseudocode 


given score. In order to fast compute an initial upper bound, Malaguti and Toth 
performed 20 runs of the greedy algorithm SEQ. 

Malaguti and Toth’s greedy approach in most case instances achieves better 
(lower) upper bounds than the greedy algorithm we took from [s^, but the 
experiments indicate that presented VNS easily decrease those (higher) upper 
bounds. Therefore, using slightly greater starting values of k could not signifi¬ 
cantly aggravate the overall optimization process. The exceptions of this ’’rule” 
can appear in some cases of small instances, where Malaguti and Toth’s greedy 
algorithm achieves nearly best known solutions, so iteration process needs to 
decrease the upper bound only for few values. In these cases, our algorithm 
needs some more time, since it starts with higher upper bounds. 

In the most recent papers ([l^, [23, 241), a constructive heuristic is not used. 
The proposed algorithms simply start with the best known values of k from the 
literature as the starting values and try to construct the feasible fc-coloring. If 
succeeds, i.e. the legal fc-coloring is found, the algorithms decrease the value 
fc by 1 and try to find the fc — 1 coloring. This iterative process stops when 
no legal coloring can be found. Although this approach appears to be very 
successful and can speed up the overall process (because the algorithms do not 
spend any time to construct some starting solution and decrease it many times 
in the optimization process), we still decided to follow the approach used in 
[iflj l and [1^: our basic approach is to construct an initial solution by greedy 
algorithm and decrease it in the iteration process, while the stopping criteria 
are not satisfied. 







3.2. The initialization and the objective function 

For the given graph G = {V,E), the solution is represented by an integer 
array of the dimension n, n = \ V\. After the upper bound (UB) is determined 
by the greedy approach, the VNS starts with the constructing the initial so¬ 
lution. Each vertex is randomly colored by a color chosen from the interval 
[1,UB — 1]. From this representation, it is obvious that the algorithm deals 
both with feasible and infeasible solutions. So, in the objective function, each 
solution is a subject of the evaluation, where the objective value is proportional 
to the ” degree of infeasibility”. In order to construct the appropriate function, 
we took into account not only the number of conflicts (pairs of adjacent vertices 
i and j for which a conflict appears) but also the ’’degree” of each conflict, i.e. 
the difference between the given distance and the distance between two assigned 
colors. Similar approach can also be seen for example in (^ [^. 

For the given solution x, represented by the array [c(l), c(2),..., c(n)], where 
c(i), i = l..n is the color assigned to the vertex i, the objective function Calcu¬ 
late () is defined as follows: 

Calculate{x) = E max{0,d(i,j) - |c(i) - c(j)|} (1) 

where d(i,j) is the given distance between the vertices i and j. From the 
equation[T]it can be seen that the infeasible solutions are penalized by increasing 
the value of the objective function, if the distance conflicts appear. In a case 
when there are no distance conflicts, the objective function is equal to 0 and in 
that case the solution is feasible i.e. the legal coloring is found. Additionally 
to the calculating the objective function, in the procedure Calculate() for each 
vertex v we remember the value conflicts{v) containing the value of distance 
conflicts of V. More precisely, this value is calculated by the formula: 

conflicts{v) = E max{0,d{v,j)-\c{v)-c{j)\} (2) 

These values, obtained for each vertex v G V, take places for the arrangement 
of the vertices, performed inside the VND procedure. 

3.3. Neighborhoods and Shaking procedure 

The shaking procedure creates a new solution x', {x' G Nk(x)), which is 
based on the current best solution x. 

In order to define the fc-th neighborhood we use the following procedure. 
Some k vertices from V are chosen randomly and for each chosen vertex, the 
color is randomly changed to some other color from the interval [l,max_color], 
where max-color is the maximal color used in the coloring. The solution x', 
formed in the described way is the subject of the further improvements in the 
VND. 

In the algorithm, the value kmin is set to 2. Experiments show that the 
algorithm achieves best performances for kmax = 20. 
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3.4- Variable neighborhood descent 

After the solution x' is obtained by the shaking procedure, a series of special 
designed local search procedures is called. In each call, the ordering of vertices, 
which are the subject of change in the VND procedure, is determined. This 
ordering enables the vertices with more conflicts to earlier become the subject 
of improvement. This number of conflicts is calculated for each vertex each time 
when the function Calculate() is called (see the subsection 13.21) . This strategy 
avoids the unnecessary steps in the local search, by increasing the probability 
that the solution will be earlier improved if the ’’worse” vertices are considered 
first, rather than if the vertices are considered sequentially without ordering. At 
the other hand, experiments show that this criterion is not enough for getting 
a quality ordering, since many vertices, after a few iterations, have the same 
number of conflicts (in the sense of the value conflicts{v)). So, additional 
criteria are involved. Firstly, we started from the fact that the colors used in 
BCP are not equal among themselves (like in the case of VCP). In fact, the 
colors which are closer to the middle of the interval [1, nc] (nc is total number 
of used colors) are harder to be replaced by another, comparing to the colors 
far off-the middle. So, if two vertices u and v have the same number of conflicts 
(conflicts{u) = con/Zicts(u)), we firstly choose u if |nc/2 — c(it)| < |nc/2 —c(u)|, 
otherwise we firstly choose v. 

In order to further improve the behaviour of the VND, we involve one ad¬ 
ditional criterion for the arrangement of the vertices, in cases when first two 
criteria do not make any difference between vertices. The third criterion is 
based on two components: for each vertex v we calculate the value weights{v), 
as the sum of the weights (distances) of the edges incident with v. Addition¬ 
ally, for the vertex v, we take into account the maximal edge weight (dis¬ 
tance) for the vertex v {maxw(v)). Finally, the third criterion is calculated 
by the geometric mean of these two values, i.e. in cases when first two cri¬ 
teria do not determine the priority, the vertex u is chosen before the vertex 
V if weights{u) ■ maxw(u) > \Jweights{v) ■ maxw{v). Otherwise, the VND 
firstly choose v. Although the third criterion influence on the ordering relatively 
rarely, the experiments show that it refines the ordering of vertices in a good 
direction and improves the overall VND. 

Beside the main experiments described in the sections 14.11 and 14.21 in order 
to justify the usage of these three criteria, additional experiments are performed. 
In the experiments (described in the section lT3l) . various combinations of the 
mentioned criteria are considered. The obtained results indicate that the ap¬ 
proach of using all three criteria generally provides better results than variants 
in which some of these criteria are omitted. 

The pseudocode of the VND procedure is shown in the figure [5] 

After the objective function for the solution x" is calculated by the function 
ObjF() and the array of vertices is arranged by using the mentioned criteria, the 
VND iteratively chooses the vertices from that array. For the selected vertex u, 
the VND is trying to find ’’better” color in the following way. If the vertex v is 
currently colored by the color c(u), we first ’’uncolor” that vertex, also removing 
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Data: x,x' 
impr <r- true; 
x" ^ x'; 
while impr do 
impr <r- false; 
objval <— ObjF(a;"); 
vertices ^ qsort(vertices,criteria); 
foreach vertex v G vertices do 
uncolor (v); 

remove_conflicts(a;", v); 
foreach color c do 
recolor(v,c); 

calculate_new_conflicts (x"); 

end 

x" = find_best_recoloring(x"); 
if (ObjF(x'f)==0) then 

//feasible coloring is found; 
x” ^ x'l; 
impr p- true; 
xs x'l; 

k* max_color(x/); 
remove Jast_color(x"); 
else 

if (ObjF(xl)< objval ) then 
//improvement happened; 
x" ^ xl; 
impr p- true; 
objval = ObjF(x/); 
end 
end 
end 
end 

Figure 2: VND pseudocode 
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the conflicts which arise by using that color for the vertex v. It should be noted 
that we do not need to calculate the objective function from the beginning, since 
only conflicts related to the chosen vertex Influence on the total sum. Therefore, 
after only these conflicts are removed, we iteratively color the vertex v by all 
other colors, trying to find better coloring. Simultaneously, we calculate the new 
conflicts, which appear by new coloring and remember these values (conflicts). 
After we tried all the colors for the vertex v, we choose the one, which generates 
the least total sum of conflicts. In that way, we get a new temporary solution 
x'l and three possibilities can happen: 

• The objective value of x'l is equal to 0: That means that the VND was 
totally successful and not just Improved the objective value, but gave the 
feasible coloring with the less number of colors than the previous best 
one. In this case, we remember that coloring, set xs = x'l and continue 
the search as follows: we set up a new current solution x" = x'l and replace 
the maximal color in x" by some other, randomly chosen color. The VND, 
l.e. the Improvement process Is then applied again on the current x" and 
the next vertex. It should be noted that if this case happen, the maximum 
number of colors may be decreased by more than one In one such step. 

• The objective value of x'l is greater than 0, but still less than the starting 
one: we hold the changes arisen in the current step (we set x" = x'l) and 
continue the VND with the next vertex. 

• The objective value of x'l Is greater or equal than the starting one: we 
continue the VND with x" and the next vertex, without any change. 

The VND finishes when no more improvements can be done and the algorithm 
analysis the results of the VND in the procedure Compare: In cases when 
number of colors used In the solution x" is less than In the solution x, or if the 
objective value of x" Is less than of the x, the currently best solution x gets the 
value x". If the objective values of the two solutions x and x" are the same, 
then X = x" is set with probability Pmove and the algorithm continues the search 
with the same neighborhood. In all other cases, the search is repeated with the 
same x and the next neighborhood. 

Although some VNS implementation allows low values oipmove, even pmove = 
0 (see for example [s^), the experiments Indicate that the proposed algorithm 
for solving BCP achieves best performances for the value Pmove = 0.5. This can 
be explained by the integer nature of the objective function and the fact that 
there is a huge number of solutions with the same objective value (especially 
when the objective value Is decreased to I). In general, higher values of Pmove 
increase the probability of the appearance of cycles. For the value Pmove = 0.5, 
the probability of the appearance of the cycle of the length I is equal to 0.5* and 
it is rather low, even for the smaller values of 1. 
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4. Experimental results 


This section presents the experimental results, which show the effectiveness 
of the proposed VNS. All the tests are carried out on the Intel 17-4770 CPU 
@3.40 GHz with 8 GB RAM and Windows 7 operating system. For each ex¬ 
ecution, only one thread/processor is used. The VNS is implemented in C 
programming language and compiled with Visual Studio 2010 compiler. 

For all experiments we used standard set of GEOM instances, which con¬ 
sists of 33 geometric graphs generated by Michael Trick, available in 13|. In 


each GEOM instance, points are generated in a 10,000 by 10,000 grid and are 
connected by an edge if they are close enough together. Edge weights are in¬ 
versely proportional to the distance between vertices, which simulates the real 
situations where closer adjacent vertices stronger interference each other. This 
set contains three types of graphs: for each dimension, one sparse (GEOMn) 
and two denser graphs (GEOMna and GEOMnb) are given. The instances are 
originally generated for BGP, but they are also transformed for solving BMCP, 
by introducing the weights of vertices. Por BMCP, vertex weights are uniformly 
randomly generated, between 1 and 10 for sets GEOMn and GEOMna, and 
between 1 and 3 for set GEOMnb. 

Since there are several state-of-the-art algorithms in the literature, we could 
not establish a unique set of algorithm parameters to make a complete fair com¬ 
parison to all these methods. Anyway, we decided to follow similar conditions 
from the most recent and most successful approaches [2S[23,|2i,|2i|. To check 
the speed of the our computer CPU, we used a standard benchmark program 
(dfmax), together with a benchmark instance (R500.5) also used in the reference 
works. For this instance, we report the computing time of 8 seconds, which is 
similar to the times reported in other recent works. Similar to the approach 
in [i^, we set the timeout limits to 2 hours for BCP instances and 4 hours for 
MBCP instances. The algorithm stops if it cannot solve the coloring within the 
time limit. For each instance, we performed 30 independent executions, which 
is also the case in 12211. 


^.1. Experimental results on BCP instances 

This section reports the experimental results obtained on the first set of 
instances, related to the BCP problem. Table [T] provide the results obtained 
by the presented VNS. The data in the table are organized as follows: first 
three columns contain the instance name, number of vertices (|U|) and number 
of edges (|i?|). The fourth column contain the best known result from the 
literature. The next five columns contain data related to the VNS: column 
(fc*) contains the best found result, the average result (column avg) obtained 
in 30 runs, the total average execution time in seconds needed to achieve the 
presented best result (column time), the hit rate {Nhn), as well as the column 
{k* — best), which contains the difference between the best result obtained by the 
VNS and the previous best-known result from the literature. A new best result 
is marked in bold. From the Table [H it can be seen that for the most of the 
smaller instances up to 60 vertices, the proposed VNS algorithm achieves best 
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Table 1: Results of the VNS obtained on BCP instances 


Instance 

|V| 

\E\ 

best 

fc”' 

avg 

time[s] 


k'^ — best 

GEOM20 

20 

40 

20 

21 

21 

0.00643 

30/30 

1 

GEOM20a 

20 

57 

20 

20 

20 

0.01507 

30/30 

0 

GEOM20b 

20 

52 

13 

13 

13 

0.0053 

30/30 

0 

GEOM30 

30 

80 

27 

28 

28 

0.01023 

30/30 

1 

GEOMSOa 

30 

111 

27 

27 

27 

0.086 

30/30 

0 

GEOM30b 

30 

111 

26 

26 

26 

0.00817 

30/30 

0 

GEOM40 

40 

118 

27 

28 

28 

0.01657 

30/30 

1 

GEOM40a 

40 

186 

37 

37 

37 

1.21 

30/30 

0 

GEOM40b 

40 

197 

33 

33 

33 

2.9218 

30/30 

0 

GEOM50 

50 

177 

28 

28 

28 

0.03407 

30/30 

0 

GEOMSOa 

50 

288 

50 

50 

50 

4.2695 

30/30 

0 

GEOMSOb 

50 

299 

35 

35 

35.0333333 

716.35247 

29/30 

0 

GEOM60 

60 

245 

33 

33 

33 

0.15027 

30/30 

0 

GEOMeOa 

60 

399 

50 

50 

50 

23.6351 

30/30 

0 

GEOMeOb 

60 

426 

41 

41 

41.8 

6430.65223 

7/30 

0 

GEOM70 

70 

337 

38 

38 

38 

0.32483 

30/30 

0 

GEOMTOa 

70 

529 

61 

61 

61.0666667 

1513.0634 

28/30 

0 

GEOMTOb 

70 

558 

47 

48 

49.2333333 

7702.62073 

2/30 

1 

GEOM80 

SO 

429 

41 

41 

41 

2.59037 

30/30 

0 

GEOMSOa 

80 

692 

63 

63 

63.3666667 

3487.41707 

21/30 

0 

GEOMSOb 

80 

743 

60 

60 

61.8666667 

7851.79623 

2/30 

0 

GEOM90 

90 

531 

46 

46 

46 

1.5592 

30/30 

0 

GEOMOOa 

90 

879 

63 

63 

64.1333333 

8082.11243 

5/30 

0 

GEOM90b 

90 

950 

69 

71 

72.7 

7627.29127 

4/30 

2 

GEOMIOO 

100 

647 

50 

50 

50 

320.8932 

30/30 

0 

GEOMlOOa 

100 

1092 

67 

68 

69.5666667 

8227.2824 

4/30 

1 

GEOMlOOb 

100 

1150 

72 

73 

75.3666667 

7945.8779 

4/30 

1 

GEOMllO 

110 

748 

50 

50 

50.0333333 

943.76663 

29/30 

0 

GEOMllOa 

110 

1317 

71 

71 

72.6333333 

9152.86497 

2/30 

0 

GEOMllOb 

110 

1366 

77 

79 

81.0666667 

8829.13587 

2/30 

2 

GEOM120 

120 

893 

59 

59 

59.0333333 

615.18443 

29/30 

0 

GEOM120a 

120 

1554 

82 

82 

84.2666667 

9112.61543 

2/30 

0 

GEOM120b 

120 

1611 

84 

86 

87.9333333 

9159.90177 

2/30 

2 


known solutions, with the exception of three instances, for which in 1^ better 
results were reported. It should be mentioned that other state-of-art algorithms 
could not achieve these results from 0 either. Since heuristic methods do 
not guarantee the optimality of the obtained solutions, it could be interesting 
if an exact method appears to check these three solutions. For the rest of the 
instances (total of 18 middle and larger instances), VNS achieves 12 best results 
and in 6 out of 18 cases VNS is not able to reach previously best-known solution 
from the literature. 

The proposed VNS solves the sparse graphs (GEOMn) easier than the dense 
ones (GEOMna and GEOMnb). From the column Nhu, it can bee seen that the 
ratio between the number of runs of the VNS when the previously best-known 
solution is reached and the total number of runs (30) is rather high for sparse 
graphs. 

Table m describes the comparison of our approach to the state-of-the-art 
methods. The data are organized as follows: The first column contain the in¬ 
stance name. The rest of the table contain data related to the approaches: the 
best result of the Discropt general heuristics (DGH) presented by Phan and 
Skiena 0 (the execution time is not reported), Prestwich’s forward checking 
colouration neighborhood search (FGNS) from [1^, Malaguti and Toth’s evolu- 
tionary algorithm (EA) from , the MITS algorithm presented by Lai and Lu 
from [^, path relinking (PR) algorithm from [ 2 ^ and learning hybrid-based 
search (LHS) from [^. All these results are extracted from [^. Last three 
columns contain best results, the execution times for the proposed VNS and the 
difference between best VNS and preyious best known result. It should be noted 
that the execution times are yery different and are not reasonable comparable. 
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since the algorithms achieve different k- colorings and the algorithms were run 
on the machines with different CPU speeds. 


Table 2: Comparison of the proposed VNS algorithm to other reference works on BCP 


instances 


Instance 

best 

DGH ^ 

PCNS 

k titles] 

EA 

[2^ 

time[s] 

MITS [2^ 
k time[s] 

PR 

k 

time[s] 

LHS 

k 

i [24] 
time[s] 

k* 

VNS 

time[s] 

diff 

GEOM20 

20 

20 

21 

0 

21 

0 

- 

- 

21 

0 

21 

0 

21 

0.00643 

1 

GEOM20a 

20 

20 

20 

0 

20 

0 

20 

0 

20 

0 

20 

0 

20 

0.01507 

0 

GEOM20b 

13 

13 

13 

0 

13 

0 

13 

0 

13 

0 

13 

0 

13 

0.0053 

0 

GEOM30 

27 

27 

28 

0 

28 

0 

- 

- 

28 

0 

28 

0 

28 

0.01023 

1 

GEOM30a 

27 

27 

27 

0 

27 

0 

27 

0 

27 

0 

27 

0 

27 

0.086 

0 

GEOM30b 

26 

26 

26 

0 

26 

0 

26 

0 

26 

0 

26 

0 

26 

0.00817 

0 

GEOM40 

27 

27 

28 

0 

28 

0 

- 

- 

28 

0 

28 

0 

28 

0.01657 

1 

GEOM40a 

37 

38 

37 

0 

37 

0 

37 

0 

37 

0 

37 

0 

37 

1.21 

0 

GEOM40b 

33 

36 

33 

0 

33 

0 

33 

0 

33 

0 

33 

0 

33 

2.9218 

0 

GEOM50 

28 

29 

28 

0 

28 

0 

28 

0 

28 

0 

28 

0 

28 

0.03407 

0 

GEOMSOa 

50 

54 

50 

2 

50 

0 

50 

0 

50 

0 

50 

0.1 

50 

4.2695 

0 

GEOM50b 

35 

40 

35 

0 

35 

0 

35 

3 

35 

1 

35 

1.2 

35 

716.35247 

0 

GEOM60 

33 

34 

33 

0 

33 

0 

33 

0 

33 

0 

33 

0 

33 

0.15027 

0 

GEOMeOa 

50 

54 

50 

1 

50 

0 

50 

1 

50 

0 

50 

0.1 

50 

23.6351 

0 

GEOMeOb 

41 

47 

43 

0 

41 

29 

41 

277 

41 

105 

41 

214.7 

41 

6430.65223 

0 

GEOM70 

38 

40 

38 

0 

38 

0 

38 

0 

38 

0 

38 

0 

38 

0.32483 

0 

GEOM70a 

61 

64 

62 

2 

61 

12 

61 

45 

61 

47 

61 

23.7 

61 

1513.0634 

0 

GEOM70b 

47 

54 

48 

1 

48 

52 

47 

8685 

47 

6678 

47 

665.4 

48 

7702.62073 

1 

GEOM80 

41 

44 

41 

0 

41 

0 

41 

0 

41 

0 

41 

0.1 

41 

2.59037 

0 

GEOMSOa 

63 

69 

63 

12 

63 

150 

63 

21 

63 

12 

63 

6.6 

63 

3487.41707 

0 

GEOMSOb 

60 

70 

61 

0 

60 

145 

60 

322 

60 

191 

60 

19.9 

60 

7851.79623 

0 

GEOM90 

46 

48 

46 

3 

46 

0 

46 

0 

46 

0 

46 

0 

46 

1.5592 

0 

GEOM90a 

63 

74 

64 

2 

63 

150 

63 

230 

63 

191 

63 

23.8 

63 

8082.11243 

0 

GEOM90b 

69 

83 

72 

2 

70 

1031 

69 

20144 

69 

23560 

69 

779.2 

71 

7627.29127 

2 

GEOMIOO 

50 

55 

50 

0 

50 

2 

50 

2 

50 

2 

50 

1 

50 

320.8932 

0 

GEOMlOOa 

67 

84 

68 

9 

68 

273 

67 

11407 

67 

5556 

67 

1557.4 

68 

8227.2824 

1 

GEOMlOOb 

71 

87 

73 

15 

73 

597 

72 

24561 

72 

41832 

71 

2038.6 

73 

7945.8779 

2 

GEOMllO 

50 

59 

50 

4 

50 

3 

50 

2 

50 

5 

50 

1.3 

50 

943.76663 

0 

GEOMllOa 

71 

88 

73 

7 

72 

171 

72 

1529 

71 

5140 

71 

2218.7 

71 

9152.86497 

0 

GEOMllOb 

77 

87 

79 

2 

78 

676 

78 

24416 

78 

18136 

77 

2598.7 

79 

8829.13587 

2 

GEOM120 

59 

67 

60 

4 

59 

0 

59 

1 

59 

2 

59 

0.5 

59 

615.18443 

0 

GEOM120a 

82 

101 

84 

4 

84 

614 

82 

34176 

82 

62876 

82 

171.1 

82 

9112.61543 

0 

GEOM120b 

84 

103 

86 

9 

84 

857 

- 

- 

85 

66301 

84 

3568.7 

86 

9159.90177 

2 


From the Table [21 it can be seen that the VNS algorithm is competitive 
with other methods. VNS achieves total of 24 previous best-known solutions. 
It is evident that the results achieved by older methods from 141 and [l^ are 
improved by the newer ones, with the exception of three small instances, for 
which the best results are reported in 141. The VNS and EA from give 
28 same best results, in 3 cases EA is better and VNS is better in two cases. 
Comparing to the MITS from for VNS and MITS 23 equal best results are 
reported, in 5 cases MITS is better, in one case VNS is better. In four cases, 
no solution is reported in [2^. VNS achieve 27 same solutions as the two most 
recent methods PR and LHS, and in 6 cases PR and LHS are better than VNS. 


4-2. Experimental results on BMCP instances 

Algorithm developed for BCP instances can also be applied for solving 
BMCP, after the implicit transformation of each BMCP to BCP instance. From 
the experimental results presented in this section we see that our VNS achieve 
many previously known best known results and two new best ones. 

Table O provides the results obtained on BMCP instances. The Table |3| is 
organized similar to the case of BCP instances: first four columns contain the 
instance name, number of vertices (|U|) and number of edges (|E|) and the best 
known result from the literature. The next five columns contain data related to 
the VNS: column (k*) contains the best found result, the average result (column 
avg) obtained in 30 runs, the total average execution time in seconds needed to 
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Table 3: Results of the VNS obtained on BMCP instances 

Instance| V|| E \best_avg_time[s ]^hit fc* — best 


GEOM20 

GEOM20a 

GEOM20b 

GEOM30 

GEOM30a 

GEOM30b 

GEOM40 

GEOM40a 

GEOM40b 

GEOM50 

GEOMSOa 

GEOMSOb 

GEOM60 

GEOM60a 

GEOM60b 

GEOM70 

GEOM70a 

GEOM70b 

GEOM80 

GEOMSOa 

GEOMSOb 

GEOM90 

GEOM90a 

GEOM90b 

GEOMIOO 

GEOMlOOa 

GEOMlOOb 

GEOMllO 

GEOMllOa 

GEOMllOb 

GEOM120 

GEOM120a 

GEOM120b 


20 40 

20 57 

20 52 

30 80 

30 111 

30 111 

40 118 

40 186 

40 197 

50 177 

50 288 

50 299 

60 245 

60 399 

60 426 

70 337 

70 529 

70 558 

80 429 

80 692 

80 743 

90 531 

90 879 

90 950 

100 647 

100 1092 

100 1150 

110 748 

110 1317 

110 1366 

120 893 

120 1554 

120 1611 


149 149 

169 169 

44 44 

160 160 

209 209 

77 77 

167 167 

213 213 

74 74 

224 224 

311 311 

83 83 

258 258 

353 353 

113 114 

266 267 

465 463 

115 116 

379 379 

357 355 

138 138 

328 329 

372 373 

142 142 

404 404 

429 429 

153 156 

375 375 

478 480 

201 202 

396 396 

536 539 

187 190 


149 

169 

44 

160 

209 

77 

167.0333333 

214.1666667 

74.0333333 

224.1 

314.2333333 
84.2 

258.4333333 

355.3 

115.7 

267.7 

465.5333333 

118.2 
381.0333333 

358.5333333 
139.0666667 
330.8333333 

374.9 

145.0333333 

405.9333 

431.9333333 

159.3666667 

376.2333333 
484 

203.4 

398.5 
545 

192.5 


54.2701 

3289.39977 

0.02153 

5.7033 

4123.7125 

1.2923 

2107.13567 

13192.2284 

1821.8568 

1671.3315 

21685.93933 

14260.65607 

5780.4842 

23989.30317 

16381.0242 

12384.2772 

20686.36647 

16783.2271 

16538.80857 

29208.93163 

13704.121 

18760.85243 

24087.21087 

19996.89127 

14816.92453 

35663.42977 

24776.3847 

22997.44417 

46954.99173 

22076.94307 

17919.32823 

59717.3869 

22835.9246 


30/30 0 

30/30 0 

30/30 0 

30/30 0 

30/30 0 

30/30 0 

29/30 0 

5/30 0 

29/30 0 

27/30 0 

2/30 0 

5/30 0 

19/30 0 

2/30 0 

2/30 1 

14/30 1 

2/30 -2 

1/30 1 

2/30 0 

1/30 -2 

9/30 0 

2/30 1 

2/30 1 

2/30 0 

6/30 0 

2/30 0 

1/30 3 

8/30 0 

1/30 2 

4/30 1 

8/30 0 

3/30 3 

3/30 3 


achieve the presented best result (column time), the hit rate (Nhu), as well as 
the column (fc* — best), which contains the difference between the best result 
obtained by the VNS and the previous best-known result from the literature. 
In column fc*, new best results are marked in bold. 

Data in Table [3] indicate that VNS achieves 21 best known results and in 
2 cases VNS obtains new best results. In 10 cases, VNS could achieve results 
close to the best ones. For all smaller BMCP instances (up to 50 vertices), VNS 
obtains all previous best-known results, with a relatively high hit rate for most 
of them. VNS achieves 2 new best results, for two middle instances GEOMfOa 
and GEOMSOa. For 9 large-sized instances, (100-120 vertices), VNS achieves 
previously known best results in 4 cases and in 5 cases VNS achieves nearly best 
results. 

TablelHshows the comparative results obtained by the state-of-the-art meth¬ 
ods from the literature and the proposed VNS. The first column of the table 
contains the instance name. The next five blocks of two columns contain best 
results and execution times of the five recent and most successful methods from 
the literature: FGNS [l^ by Prestwich, Malaguti and Toth’s EA from (20l |. 
the MITS algorithm from [22| presented by Lai and Lu, path relinking (PR) 
algorithm from [23| and learning hybrid-based search (LHS) from [ 2 ^. Like 
in the case of BCP, all these results are also extracted from [^. Last three 
columns contain best results, the execution times for the proposed VNS and the 
difference between best VNS and previous best-known results. 

Prom the Table|4]it can be seen that FCNS, EA, MITS and PR could achieve 
most of the best known solutions for small instances up to 50 vertices, and LHS 
and VNS achieve all of the best known solutions for small instances. MITS and 
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Table 4: Comparison of the proposed VNS algorithm to other reference works on BMCP 


instances 


Instance 

GEOM20 

GEOM20a 

GEOM20b 

GEOM30 

GEOM30a 

GEOM30b 

GEOM40 

GEOM40a 

GEOM40b 

GEOM50 

GEOMSOa 

GEOM50b 

GEOM60 

GEOM60a 

GEOMGOb 

GEOM70 

GEOM70a 

GEOM70b 

GEOM80 

GEOMSOa 

GEOMSOb 

GEOM90 

GEOM90a 

GEOM90b 

GEOMIOO 

GEOMlOOa 

GEOMlOOb 

GEOMllO 

GEOMllOa 

GEOMllOb 

GEOM120 

GEOM120a 

GEOM120b 


best k time k time 


149 149 

169 170 

44 44 

160 160 

209 214 

77 77 

167 167 

213 217 

74 74 

224 224 

311 323 

83 86 

258 258 

353 373 

113 116 

266 277 

465 482 

115 119 

379 398 

357 380 

138 141 

328 339 

372 382 

142 147 

404 424 

429 461 

153 159 

375 392 

478 500 

201 208 

396 417 

536 565 

187 196 


4 149 18 

2 169 9 

0 44 5 

0 160 1 

11 210 954 

0 77 0 

1 167 20 

299 214 393 

4 74 1 

1 224 1197 

51 316 4675 

1 83 197 

77 258 139 

10 357 8706 

12 115 460 

641 272 1413 

315 473 988 

55 117 897 

361 388 132 

109 363 8583 

37 141 1856 

44 332 4160 

13 382 5334 

303 144 1750 

7 410 3283 

26 444 12526 

367 156 3699 

43 383 2344 

29 490 2318 

5 206 480 

9 396 2867 

41 559 3873 

3 191 3292 



149 2 

169 15 

44 0 

160 0 

209 10 

77 0 

167 0 

213 328 

74 2 

224 8 

314 40373 

83 1200 

258 19 

356 38570 

113 104711 

270 7602 

467 38759 

116 213545 

381 212213 

361 41235 

139 255 

330 4022 

375 10427 

144 211366 

404 40121 

442 381 

156 213949 

381 183 

488 926 

204 944 

554 1018 

189 213989 



149 1 

169 7 

44 0 

160 0 

209 26 

77 0 

167 1 

213 133 

74 4 

224 2 

312 270860 

83 723 

258 23 

354 34580 

113 63579 

266 130844 

466 6952 

116 26110 

380 34493 

358 41772 

138 705 

328 134941 

372 282456 

144 14648 

404 16355 

436 9108 

156 86308 

375 25401 

482 9819 

201 47653 

396 15341 

539 45147 

189 14371 



149 1.8 

169 0.5 

44 0 

160 0.1 

209 16.2 

77 0 

167 0.2 

213 9 

74 1.5 

224 0.3 

311 1452.6 

83 72.1 

258 1.3 

353 9007.1 

113 910.7 

266 2534 

465 36604.9 

115 3640.7 

379 357.8 

357 43403 

138 46.5 

328 162.2 

372 16782.1 

142 7680.8 

404 64.9 

429 78363.1 

153 10840.1 

375 1598.8 

478 49457.1 

201 5388.4 

396 626.1 

536 69518.6 

187 8025.8 


k* 

149 

169 

44 

160 

209 

77 

167 

213 

74 

224 

311 

83 

258 

353 

114 

267 

463 

116 

379 

355 

138 

329 

373 

142 

404 

429 

156 

375 

480 

202 

396 

539 

190 


VNS 

time[s] diff 
54.2701 0 

3289.39977 0 

0.02153 0 

5.7033 0 

4123.7125 0 

1.2923 0 

2107.13567 0 

13192.2284 0 

1821.8568 0 

1671.3315 0 

21685.93933 0 

14260.65607 0 

5780.4842 0 

23989.30317 0 

16381.0242 1 

12384.2772 1 

20686.36647 -2 

16783.2271 1 

16538.80857 0 

29208.93163 -2 

13704.121 0 

18760.85243 1 

24087.21087 1 

19996.89127 0 

14816.92453 0 

35663.42977 0 

24776.3847 3 

22997.44417 0 

46954.99173 2 

22076.94307 1 

17919.32823 0 

59717.3869 3 

22835.9246 3 


PR fails to achieve best known solution only for one small instance (GEOM50a). 

For the middle-size instances (60 to 90 vertices), the situation is slightly 
different; FCNS and FA could achieve only one best known solution (instance 
GFOM60), while MITS could achieve two best known solutions for two instances 
with 60 vertices. For other middle-size instances, MITS achieves nearly best 
known solutions. For this class of instances, PR could achieve 6 best known 
solutions and in other cases is close to the best ones. LHS achieves all best 
known solutions except in ten cases and in two cases close to the best ones. 
VNS achieves five best known solutions, in five cases VNS achieve solutions 
close to the best ones and in two cases VNS finds new best solutions. 

For the large-size instances, PR could not find any best known solution, while 
FA and MITS achieve only one best known solution. PR achieves four best 
known solutions and for other large-size instances achieve solutions relatively 
close to the best ones. While LHS achieves all best knows solutions for large 
instances, VNS achieve best ones for six instances and for the rest nearly best 
solutions. 

Regarding execution times, it is obvious that they are quite different, since 
different approaches provide results with different k values. Additionally, some 
of the algorithms (like FA and VNS) spend additional time for the construction 
phase, while other algorithms (MITS, PR and LHS) start with the previous best 
known solutions. 
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4-3. Justification of the usage of the criteria in the VND procedure 

In order to further analyse the behaviour of the proposed VND procedure 
and to justify the usage of the criteria for ordering vertices, additional exper¬ 
iments are performed with three large BMCP instances (namely GEOMllO, 
GEOMllOa and GEOMllOb). Recall that in the VND procedure three criteria 
are combined for determining the ordering: the first criterium is the total num¬ 
ber of conflicts for the vertex, the second is the distance of the color assigned 
to the vertex from the middle and the third one is the geometric mean of the 
total sum of the weights of the edges incident with the vertex and the maximal 
edge distance for that vertex. In the experiment, total of eight possible differ¬ 
ent variants of the overall criterium are analysed: each of the mentioned three 
criteria is or is not used. 

Table [5] provides the experimental results obtained in all of these eight cases. 
Each case is denoted as XYZ, where X, Y and Z belongs to {0,1}, indicating 
if the criterium is used (value 1) or not used (value 0). For each case, best 
found and average results are shown. The time limit for this experiment is 
set to 3000s. From the table [SJ it can be seen that the variant proposed in 


Table 5: Comparison of the various combinations of the criteria used in the VND 



000 

001 

010 

Oil 

instance 

best 

avg 

best 

avg 

best 

avg 

best 

avg 

GeomllO 

376 

379.16667 

375 

378.53333 

375 

378.96667 

375 

378.46667 

GEOMllOa 

486 

491.23333 

487 

491.23333 

483 

487.86667 

483 

488.33333 

GEOMllOb 

201 

204.06667 

201 

204.3 

202 

204.1 

202 

203.96667 


100 

101 

110 

111 

instance 

best 

avg 

best 

avg 

best 

avg 

best 

avg 

GeomllO 

376 

378.2 

376 

378.03333 

376 

378.56667 

375 

377.9 

GEOMllOa 

482 

489.83333 

484 

488.96667 

483 

487.73333 

481 

486.8333 

GEOMllOb 

203 

204.6 

202 

204.6 

202 

204.16667 

202 

204.0333 


this paper (column III - all three criteria are used) tends to be the best one, 
although some other variants also gives good results. For the first instance 
(GEOMllO), in the variants denoted as 001, 010 and 011 the best result is also 
achieved, but the average result is worse than in the proposed variant (HI). For 
the second instance (GEOMllOa), in the proposed variant (111) both the best 
known and the average values are better than in the other cases. For the last 
instance (GEOMllOb), the first two variants achieves better best result than 
the proposed one, but the average value obtained by the proposed variant is 
better. 

These results indicate that the usage of the proposed combination of the 
criteria is justified, having in mind that some other variants can also reach high 
quality results. 


5. Conclusions 

In this paper we present the VNS algorithm for solving two generalizations of 
the vertex coloring problems: bandwidth coloring problem and multiple band¬ 
width coloring problem. Since BCP and MBGP enjoy many applications, pre¬ 
sented algorithm and achieved results are of a great interest for both theory and 
practice. 
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In the shaking procedure, an increasing number of vertices are permuted, 
forming the new solution which is subject of the further improvement in the 
VND. In the VND procedure, vertices are arranged in a way to increase the 
probability of successful recoloring. This approach of the arrangement of the 
vertices splits the local search in a series of disjoint procedures, enabling better 
choices of the vertices which are addressed to re-color. The overall criterion 
for the ordering the vertices is based on the number of conflicts of each vertex, 
combined with two additional criteria. 

The algorithm is tested on the common used instances. In the case of BCP, 
VNS achieves many previous best-known results and in the case of BMCP, 
VNS obtains two new best solutions in a reasonable time. Experimental results 
indicate that the proposed VNS becomes one of state-of-the-art methods for 
solving BCP and BMCP. 

The further investigation of this problem can include parallelization and 
running on some powerful multiprocessor system, as well as the application of 
the proposed method for solving some similar real-life problems. 
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